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The NP-complete problem of the travelling salesman (TSP) is considered in the framework of 
quantum adiabatic computation (QAC). We first derive a remarkable lower bound for the computa- 
tion time for adiabatic algorithms in general as a function of the energy involved in the computation. 
Energy, and not just time and space, must thus be considered in the evaluation of algorithm com- 
plexity, in perfect accordance with the understanding that all computation is physical. We then 
propose, with oracular Hamiltonians, new quantum adiabatic algorithms of which not only the lower 
", bound in time but also the energy requirement do not increase exponentially in the size of the input. 

■ Such an improvement in both time and energy complexity, as compared to all other existing algo- 

' rithms for TSP, is apparently due to quantum entanglement. We also appeal to the general theory of 

, Diophantine equations in a speculation on physical implementation of those oracular Hamiltonians. 

S3 

S: 

^ ■ A lower bound on the adiabatic computation time 

Right from the earlier days of quantum computation, NP-complete problems have been the sub- 
ject of investigation within standard quantum mechanics Q] as well as in non-standard theory Q . 
Then,with the advent of QAC Q they are further pursued in the form of satisfiability prob- 
lems 0,H,@ mostly with inconclusive results 0,01 (see also QAC has also been applied 



l applied 

M0, 



to the search of an unstructured database with the same time complexity as with Grover'sL 
or with a much shorter time but at a more substantial cost for the energy involved |l3l Il4lll5j| . 

QAC starts with the readily constructible ground state \gi) of an initial Hamiltonian Hj which 
is then adiabatically extrapolated to the final Hamiltonian Hp so that the ground state \g) of the 
latter which contains the information of the solution of the problem we want to solve could be 
obtained with reasonably high probability. The interpolation between Hi and Hp is facilitated 
f> \\ by a time-dependent Hamiltonian in the time interval < t < T, 



■ H(t) = f{t)H,+g{t)H P , (1) 

a _ 

either in a temporally linear manner (in which case, f(t) = (l — ~) and g(t) = or otherwise 
(but with /(0) = 1 = g(T) and f{T) = = ,g(0)), and with or without the condition f(t)+g(t) = 1 
for all t. Such a time evolution is captured by the Schrodingcr equation 

•i-H . 

idt\m) - nt) w)h (2) 
\m) - !<?/>• 

We next consider another state \<fi{t)} which satisfies a closely related Schrodinger equation 

id t \<f>(t)) = (f{t)H I+ pg(t)\) |0(t)), (3) 
10(0)) = Iff/), 

with arbitrary f3, resulting in a phase ambiguity. Clearly up to a phase, 

\m=e m \9i)- (4) 
From the difference between Eqs. and 

d t (|V(i)) - \m) = -iH(t) (\m) - W))) + ig(t)(H P - /?)!#*)), (5) 
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we have 

d t \\\m)~\m\\ 2 = -2M({m\-m\)d t (\m)-\m), 

= 2 5 (t)9f (m)\-(<t>(t)\)(H P -0)\<j>(t)}, 

< 2g(t)\\(H P -0)\m\\x\\\m)-W))l 
9t\\W))-W))\\ < g(t)\\(H P -P)\9i)l (6) 

where the first inequality is a result of the Schwarz inequality; and the second inequality from Q . 
Integrating the time variable on both sides of the last expression from to T yields 



\\(H P - (3)\9i)\\ 

The numerator of the left hand side is the distance between the final state of \ijj(T)) and 
essentially its initial state at t = 0, 

< 2. (8) 

Now, with the help of the following result 



\\(Hp-p)\ gi ) || = V(rM(H P -(3r\ gi ), 



[g I \H 2 p \g I )-2f3(g I \Hp\g I }+02, 
= J((gi\H P \ gi ) - ( 9l \Hp\ 9l )*) + OS - {gi\H P \ gi )f, 



> yJ{m\H%\gi)-{!ii\Hp\g I )* = A I E, (9) 

we have also 

wmn - Mm < 2 ^ (io) 



||(tfp-/3)|ffi)|| " min^lK^p-^lsj)!) " AjS 

This represents an upper bound whereby there is no restriction on how large the distance from 
the final state, \ip(T)), to the initial state, \4>(T)) = e^ (T) |g/), must be. 

Now, as the right hand side of Q is a monotonically increasing function of T (thanks to the 
positivity of g(r)), we can choose a minimum time T m i„ in such a way that this right hand side is 
larger than the right hand side of (|10fl . which is independent of T, 

< / ml ".9(r)dr. (11) 



AiE 



The meaning of this time T m ; n must then be derived from the meaning of the right hand side 
of I|1U|) . Since this latter quantity represents an unhindered and maximally allowable exploration 
by the state vector \ip(T)) of the whole Hilbert space, we thus see that T min is accordingly the 
minimum time for such an unhindered exploration, in order to realise in general the full potential 
of the adiabatic computation. That is, T m j n is a lower bound on the evolution time whereby the 
whole Hilbert space can be explored unhindered, so that the final ground state, in particular, may 
be obtained with some high probability. 

Of course, the computation could be run for a time T less than this T m ; n , but then sufficient 
computation time has not in general been given for a full exploration of the underlying Hilbert 
space. For any run time less than T m i n , the state |^(T)) can explore only some smaller portion of 
the Hilbert space. (As can be seen from <|10[) . this portion is defined by a restrictive constraint on 
the distance from the state at time T, \ip(T)), to its initial state, \gi) = e _I ^ T ) \4>(T)).) 

The last inequality (|llf> could also be expressed in the form of a time-energy relation, 

2 < g(0) T min AjE, (12) 

with some 6 in the range < 9 < T m ; n , thanks to the mean value theorem. This is our first main 
result El. 
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This lower bound on the computation time incorporates the initial ground state \gi) and the 
spectrum of the final Hamiltonian together in AjE, which is defined in as the energy spread 
of the initial state in terms of the final energy. The manner of the time extrapolation is further 
reflected in g(0). The condition (|12fl states that the more the spread of the initial state in energy 
with respect to the final Hamiltonian, the less the lower bound on the running time. 

It thus also emphasises the fact that energy must be considered in the running of QAC, and 
perhaps in all physical computation if they are indeed all physical in the end. It is interesting to 
further note that recent results obtained in [tj reflect well in the condition above. 

This condition (|12fl must also contain in it the information about the energy gap between the 
instantaneous ground state and the first excited state, which in turn determines a lower bound on 
the computation time as indicated by the quantum adiabatic theorem. The dependence on the 
gap is not manifest but once again hidden in A/_E, the energy spread of the initial ground state in 
terms of the energy eigenstates of the final Hamiltonian. This is not so surprising because (|T^|) is 
derived directly from the Schrodinger equation, which is also the starting point for the derivation 
of any version of the quantum adiabatic theorem. We will see in the quantum adiabatic algorithms 
of the next Section that such an implicit gap dependence is born out in the perfect agreement of 
the lower bound i|12|) with all other lower bounds obtained directly from the energy gaps, if and 
when such gaps can be obtained. 

Energy in the consideration of algorithm complexity 

We now apply our lower bound to various quantum adiabatic algorithms in the literature. 

In order to find the item labeled by |m) in an unstructured database with N items, one 
method ^3 is to employ the initial Hamiltonian Hi = 1 — \gi)(gi\, with the initial ground state 
\gi) that has equally distritbuted probability amplitudes among the N items, \gj) = ^= X^fcLi N)i 
together with the final Hamiltonian Hp = 1 — \m)(m\, which admits |m) as the non-degenerate 
ground state. In this case, AjE = ^JN — 1/N. 

With linear time extrapolation Q also employed in the adiabatic computation Q , we substitute 
gir) = ^ into (|12f> to obtain 

0(VN) < T min (13) 

This lower bound can indeed be attained in a local adiabatic search 1 1 l2f | . Its time complexity, even 
though better than a straightforward classical search, is the same as Grover's quantum search. 

Alternatively, for example, the same Hamiltonians but with different extrapolation as proposed 
in MM, 

9{r) = ^ + Vn^ (l - J) => J g(r)dr ~Tx 0(Vn), (14) 

would lead to 

0(1) < Tmin, (15) 

that is, a constant lower bound on the computation time, irrespective of the database size! This, 
however, could only be obtained at the expense of the energy being scaled as 0(\HSf) (contained 
in the term g(t) hidden in <|14|) h even if onl y fo r the intermediate time. 

Even though being known previously [13, Ej, [Hj , we want to reiterate here that energy consid- 
eration emerges naturally from QAC, as predominantly displayed in the relation l|12|) . and wish to 
re-emphasise here the inevitable fact that, because all computation is physical, energy must be an 
independent dimension, besides those of time and space, in the evaluation of algorithm complexity 

An adiabatic algorithm for TSP with non-polynomial energy resources 

We will deal with the version of TSP 0] with M cities to be visited as a search problem for 
the shortest tour among Ml possible tours, each of which connects all the cities and visits each 
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city once only. The distances between pairs of cities can be presented in a distance matrix, which 
may or may not be symmetrical, having elements d^ being the (oriented) distance from the i-th 
city to the j-th city, and having vanishing diagonal elements, da = 0. 

A tour is labeled by an integer k, k = 1, . . . , Ml, which is the rank of a permutation of (1, ... , M). 
The length of a tour is labeled by Ik , which is the sum of appropriate dy . We also denote by i max 
a constant (not much) greater than the maximum tour length. 

We state here an assumption (needed for this Section only) about the variance Em in Ik, 



Em = 



\ 



1 M! 2 

~M\ (' fe — decreases, if at all, less than exponentially in M, (16) 

' k=l 



for asymptotically large M. That is, we expect that when we add one more city to the list to be 
visited, even though the number of possibilities now shoots up to (M + 1)!, the difference/spread 
between the longest and the shortest tour lengths should not decrease exponentially. Indeed, with 
randomly distributed dij we can, thanks to the central limit theorem, show that Em ~ 0(\fM) 
for M — > oo, namely, as Em does increase, the assumption 11161 is satisfied. We suspect that 
the assumption can be proven in the Euclidean version of the TSP, which is also an NP-complete 
problem [TtJ. In fact, it would be sufficient for our purpose if the TSP with the condition l|16f) is 
provably NP-complete still. 

To employ QAC, we introduce an oracular Hamiltonian Ho, of infinite dimensions if necessary, 
which is diagonalised in the number states \n), 

, i„ i i / Wl: for < n < Ml - 1, . . 

(n\H \n) = | Zm + xj forM! < n (17) 

Possible physical implementation of this oracular/hypothetical Hamiltonian will be discussed later. 

We now need to find the ground state |no) of Ho in order to obtain the tour information from 
uq and the least tour length from the ground state energy, assumin g n o d egen eracy. Following a 
quantum adiabatic algorithm introduced for Hilbert's tenth problem [I&ll&l2(ll2Ul2l3 |. we employ 
a linear time interpolation from the following initial Hamiltonian, with annihilation (creation) 
operators a (a^), 

Hi = (a) - a*){a- a), (18) 

|q| 2 „ 

which admits as the ground state the coherent state |a) = c ~ Y^=o "^7 \ n )- The choice of the 
complex number a will be crucial for the complexity of our algorithm. 

The initial ground state |a) has a Poissonian probability distribution among the number states 
\n), with equal average occupation number and variance, n = a 2 = \a\ 2 . As the oracular Hamil- 
tonian Ho <|17fl is only non-constant for n < Ml — 1, the choice 

H 2 = 0(Af!) (19) 

would yield, in the initial state, a variance a 2 in n that covers all the Ml possibilities, and thus 
would lead to AiE in the inequality l|12|) taking the same magnitude as the measure Em of the 
spread in the tour lengths among the Ml tours, 

A/£~0(E M ). (20) 

By the assumption l|16|) and according to l|12|) . we now have a lower bound for the computation 
time which is now of the same order as the inverse of Em and thus cannot grow exponentially 
in M. That is, we would gain in the computation time complexity if such a lower bound in time 
could be achieved for the QAC. 

However, the price we must pay for such a sub-exponential, or even constant, growth is the 
exponential cost of the energy that must be supplied, as reflected in the choice (|19fl . \a\ 2 = 
O(Ml), for the initial Hamiltonian (|18f) . For much larger values of \a\ 2 , the variance a 2 would be 
concentrating in large values of n, and accordingly AiE would be vanishingly small as is evident 
from the constant eigenvalues of Ho in (|17fl for Ml < n. For much smaller value of \a\ 2 , the 
variance AiE could be concentrating in an inappropriate domain of l n which in general cannot 
lead to a sub-exponential or polynomial lower bound for the time complexity. 

The situation is then similar to other algorithms of a previous Section which gain in the time 
complexity at the expense of the energy complexity. Fortunately, we could improve upon the 
situation to reduce the energy consumption. 
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Another algorithm for TSP with polynomial energy resources 



Instead of encoding a tour by a single number k, we now employ the M-tuple (mi, . . . , 771m}- 
Each m, is a natural number and its first M values (for i ranges from to M — 1) represent the M 
cities. For these, there are M M tuples, including those with m, = rrij, which do not correspond to 
a TSP tour. In order to put those that are not a TSP tour out of reach, we introduce a counterpart 
l s of the tour length 

j _ ( tour length l s , if s corresponds to a TSP tour, , * 

s \ d 2 + l max , otherwise, ^ ' 

where we could use, for example, the following correspondence between s and the M-tuplcs 

s = 1 + mi + m 2 M + m 3 M 2 + ... + m M M {M ~ 1] . (22) 

In (|21fl . d is a random number drawn from some arbitrarily chosen distribution that has a finite 
variance independent of M, ad = C, so that for those s that are not a TSP tour, l s > i max - 
Alternatively, one could obtain the same result by choosing d 2 = for even s, and d 2 = 2i max for 
odd s (that is, d 2 = (1 - (-l) s )Z ma x)- 

We now start our QAC with the initial coherent state \tp(0)) = ®fLi\ a i), which is the ground 
state of the initial Hamiltonian 



H, 



M 

E 



(a\ - a* J (a, - on) , (23) 



which is in turn extrapolated linearly in time to an oracular Hamiltonian that is diagonalised in 
the number states \mi . . . mu) according to 

/ i£t i \ f h, for 1< s < M M , , n ,s 

(mi... m M \Ho\mi . . . m M ) = | ^ for ^ (24) 

The aim of our computation is once again to find the ground state of Hp, with the ground state 
energy being the shortest tour length. 

If we choose all the complex numbers a^s such that 

|«,| 2 = 0(M), (25) 

then the energy requirement for l|23|l is only of 0{M 2 ) 1 which is significantly lower than that of 
O(Ml) of the last Section. Furthermore, with this choice, the variance/spread in s in the initial 
state, as can be computed easily from H22fl. will be of 0(M M ). And with this spread the energy 
variance JjJJl of the initial state in terms of the energy of (|24|l will also be of the same magnitude as 
the variance in the "tour lengths" l s . That is, for large M. AjE — > 0(<Jd), which is a predetermined 
constant independent of M, in an exponentially manner because of the exponential decrease in 

the proportion of TSP tours, (M\/M M ) AI ^° c~ M /VM. 

Applying our inequality (|12f> here once again, we are led a lower bound in computation time 
approaching a constant value independent of the number of cities M. But this time the energy 
cost is only quadratic in M . This is our second main result. 

Classical simulations of this QAC is rather expensive because of the M M possibilities, instead 
of Ml as in the usual formulation of TSP. Hence, the extraordinary gain above would only be 
possible with quantum computation. We suspect that quantum entanglement, which has not been 
fully exploited in QAC up until now, is responsible for this remarkable improvement. 

The above algorithm can also be captured in the following QAC in a finite-dimensional Hilbert 
space with linear time interpolation between 

M-l M-l 

Hi = 1- \gi)(gi\, and H P = ^ ... ^ l s \mi . . . m M )(mi . . . m M \, (26) 
where \ip(0)) = \gi) = 
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Note that similar kind of oracular Hamiltonians is already employed in Grover's search [12J and 
involves some kind of non-local interactions, which may make their implementation more difficult 
than that of quantum circuits (which only require local interactions) but not impossible. However, 
the infinite-dimension formulation above (in l\'2'.>li and l]24p) could exhibit explicitly the requirement 
of energy through the a^s in (123), and may be implemented physically with quantum optical 
means, among others. 

Physical Hamiltonians from the oracular Hamiltonians via Diophantine equations 

For completeness, we outline here some ways forward, in principle, to a physical implementation 
of the oracular/hypothetical Hamiltonians l|17l I24fl . The oracular Hamiltonian (|17|) . for example, 
can be regarded as a computer program, with n being an input upon which y n — (n\Ho\n) is 
the output (without loss of generality, we choose y also to be natural numbers). According to 
the general theory of Diophantine equations |23l |24J , corresponding to such a computer program 
(Turing machine) there is a Diophantine equation which only has integral solution in (xi, . . . , Xk )> 

D(n,y n \x u ...,x K ) = Q> (27) 

for some K , if and only if y n is indeed the output for the input n. This is a manifestation of the 
computational universality of Diophantine equations. 

Following [Tsl . ITfll l2flL I2H l22l | , we could then try to implement the physical Hamiltonian 

Hp = / max (d (ala n , ala y ; a[ai, . . . , a^ K a K ^ + a\a y . (28) 

Thanks to the fact that Diophantine polynomials have only integral values, it can be shown that 
the ground state energy of this Hamiltonian is the shortest tour length we are after. This is 
the most general but certainly not the most efficient way to consider physical implementation 
of the oracular Hamiltonians. The appearance of extra K modes in ()28[) would require further 
consideration to ensure that we still have non-exponential complexity in both time and energy. 
These are outside the scope of the present paper. 

Summary and concluding remarks 

We derived a remarkable lower bound on the computation time for QAC in general as a function 
of the energy necessary for the computation. Energy, being a physical quantity and an inseparable 
component of any physical process, emerges naturally in QAC, as predominantly displayed in the 
relation (|12|l . and is thus an inevitable and a natural dimension in the complexity evaluation of 
an algorithm, in accordance with the view that all computation is physical. We then proposed 
some new quantum adiabatic algorithms, with oracular Hamiltonians, for the NP-complcte TSP. 
Both the lower bound of the computation time and the energy required of our proposal do not 
scale exponentially with the number of cities involved. Enlisting the help of general theory of 
Diophantine equations, we speculated on a general way in principle to physically implement the 
oracular Hamiltonians of the algorithms. This, however, and whether such a favourable lower 
bound in the computation time could be achieved would require further detailed investigations 
elsewhere. 
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